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Abstract 



PLh 

• ■ This paper is the continuation of our earlier paper [5], where we proved 

^~~{ \ i^' ^-order of current fluctuations across the characteristics in a class of one 

C^ ■ dimensional interacting systems with one conserved quantity. We also 

claimed two models with concave hydrodynamic flux which satisfled the 

assumptions which made our proof work. In the present note we show that 

the totally asymmetric exponential bricklayers process also satisfles these 

Cn ' assumptions. Hence this is the flrst example with convex hydrodynamics 

^ , of a model with f^' '^-order current fluctuations across the characteristics. 

^^ ■ As such, it further supports the idea of universality regarding this scaling. 

C^^ ■ Keywords: Interacting particle systems, universal fluctuation bounds, t^/'^- 

\l^ , scaling, second class particle, convexity, bricklayers process 
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1 Introduction 

- 1 — I ■ 

^^ \ It is conjectured that particle current through the characteristics of one dimen- 

sional stochastic interacting systems with one conserved quantity and concave 
or convex hydrodynamics show i^/'^-order fluctuations and Tracy- Widom type 
limit distributions in this order. Our earlier paper ^ provides a robust argu- 
ment that proves this order of the fluctuations. We refer to that paper for the 
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general framework and other results of the field. Very briefly, [5] works if one 
proves the following properties of a model (see the exact formulation therein) : 

1. a strict domination of a second class particle of a denser system on one of 
a sparser system, 

2. a non-strict, but tight, domination of a second class particle on a system 
of second class particles that are defined between the system in question 
and another system with a different density, 

3. strictly concave or convex, in the second derivative sense, hydrodynamic 
flux function of the hyperbolic conservation law obtained by the Eulerian 
limiting procedure, 

4. a tail bound of a second class particle in a(n essentially) stationary process. 

Properties [T] and [5] form what we call the microscopic concavity or convexity 
property. Arguments in [S] are worked out for the concave setting, but every- 
thing works word-for-word in the convex case. 

We try to very briefly indicate here how the argument [5] proceeds from these 
assumptions. We consider two coupled processes with different densities, and a 
density of second class particles between them. We tag one of these second class 
particles, let us denote its position by X. We also consider the position Q of 
a single second class particle that evolves on one of the two coupled processes. 
Microscopic concavity allows us a (non strict, but tight) comparison between 
the positions X and Q. Briefly, the argument now goes as follows. A deviation 
event of Q implies, via the comparison, that of X. Next we relate the deviation 
of X to a surface growth deviation by simply noticing that current of second 
class particles is just difference in surface growth. The probability of the surface 
growth deviation is controlled by the variance. An exact connection [8^ between 
the variance of the surface growth and the flrst centered moment of Q then 
"closes the loop" , and now we have a deviation bound of Q in terms of its 
first centered moment. If we close one of the two densities to the other in the 
right pace as a function of time, then constants appear in a way that allows 
to conclude the time^' '^ scaling; this is the point where strict convexity of the 
hydrodynamic fiux function is essential. 

Two examples are claimed in [5]: the asymmetric simple exclusion process 
and a totally asymmetric zero range process with jump rates that increase with 
exponentially decaying slope. In this note we prove the above properties and 
hence the t^''^ scaling for yet another system, the totally asymmetric exponential 
bricklayers process (TAEBLP). This model was introduced in [1], and its nor- 
mal fluctuations off-characteristics were demonstrated in P] (in case of general 
convex jump rates, not only exponential). 

General convex increasing rates of a totally asymmetric bricklayers process 
allow couplings that prove properties[T]and[3]above. The exponential jump rates 
have a strong enough convexity property that will allow us to show property [2l 
We do this by repeating an argument somewhat similar to the one applied to 
the concave zero range process in [S]. The idea resembles much to the concave 



case, but we include this convex case in full details (rather than listing all the 
differences from previous work) due to the complexity of the method. 

Finally, property|3]is highly nontrivial when the jump rates have unbounded 
increments. We use a coupling based on property [1] and a recent result [1] that 
asserts that a second class particle of the exponential bricklayers process per- 
forms a simple (drifted) random walk under appropriate shock initial conditions. 
It is worth noting that exponential jump rates were also of fundamental impor- 
tance in [3], this technical point being the main reason for considering this 
particular family of jump rates in this note. Indeed, this is the only point that 
prevents us from proving the result for e.g. the totally asymmetric zero range 
process with convex exponential jump rates. 

We emphasize at this point that we only consider nearest neighbor models. 
We believe that, as far as the hydrodynamic flux is strictly convex in the second 
derivative sense, the t^^^ scaling should hold for a wide class of non nearest 
neighbor dynamics as well. However, as intricate couplings and orderings of 
second class particles play a crucial role in the methods, we do not see an easy 
way to deal with the non-nearest neighbor case. 

Let us also have a few comments on explicit product invariant stationary 
distributions. In [S] we explicitly use them, as they make the arguments easier. 
The crucial points of the method are properties 2 (microscopic convexity) and 4 
(tail bound of the second class particle) above. These depend on the details of 
the models, and the few known examples for which they could be proved indeed 
have product stationary distributions. Therefore we have not investigated how 
the arguments in [5] and the present note could be generalized to the case of 
other types of stationary distributions. We again believe that once microscopic 
convexity and the tail bound were proved, the remainder of the argument could 
be generalized and the scaling would remain valid for many models with non 
product stationary distributions as well. 

The case of exponential jump rates was constructed in [6]. The results of 
the note [5] are used by |S] . Those require strong construction results which are 
not provided by fF; and therefore, to our knowledge, are not available. To close 
that gap, we reproduce the results of [5] here for the TAEBLP. 

The organization of this paper is the following: we repeat the introduction 
of the model, the fluctuation results and conclusions of [5] , and the definition of 
the microscopic convexity property in Section [2] We construct the four process 
coupling and prove the microscopic convexity property in Section [3l Finally we 
show how to use ^ to prove property |4] in Section |4l The result of the note [8] 
is reproduced in the Appendix. 

2 The model, properties and results 

2.1 The model 

The model we discuss is the totally asymmetric exponential bricklayers process 
(TAEBLP) introduced in [I], and also treated in [3] and gj. The model is a 



member of the class in [S] , here is a brief definition. The process describes the 
growth of a surface which we imagine as the top of a wall formed by columns 
of bricks over the interval (i, i + 1) for each pair of neighboring sites i and 
i + 1 of Z. The height hi of this column is integer-valued. The components 
of a configuration w G fi are the negative discrete gradients of the heights: 
uji — hi-i — hi e Z. The configuration space is therefore 

At times it will be convenient to have notation for the increment configuration 
5_i £ U, with exactly one nonzero entry equal to 1: 

{1, for i = 7, 
f -J- 
0, for I ^ J. 

Bricklayers processes are characterized by a function / : Z — > M+. We 
only consider the totally asymmetric nearest neighbor case here, in which only 
deposition of bricks in the following way are allowed: 

/oo\ (w^w.+i)— ^ (w, -1, Wi+i + l)l , fi ^ , fi ^ 

(2.2) , , -, ) with I'ate /(w^j + /(-cj^+i). 

hi — > hi + 1 J 

Conditionally on the present state, these moves happen independently at all 
sites i. We can summarize this information in the formal infinitesimal generator 
L of the process w(-): 

iei, 

A jump results in the new configuration uj}"^' '+i' defined by 

jj for j y^ i, i + 1, 

^^(^,^+l)y^ = <( a;, -1 for j = i, 

Wj+i + 1 for j = i + 1. 

L acts on bounded cylinder functions (^ : f2 — > M (this means that ip depends 
only on finitely many w^-values). The additive form of the rates gives rise to 
the bricklayers representation: at each site i stands a bricklayer who places a 
brick on the column on his left with rate f(—uJi) and on the one on his right 
with rate f{uji). 

Thus we have a Markov process {w(t) : t £ M+} of an evolving increment 
configuration and a Markov process {h{t) : t e K+} of an evolving height 
configuration. The initial increments u;(0) specify the initial height h{0) up to 
a vertical translation. We shall always normalize the height process so that 
ho{0) = 0. 

Attractivity of the process is essential for this paper. This is achieved by 
assuming that / is nondecreasing. 



Finally, stationary translation-invariant product distributions for a;(-) are 
ensured by /(z) • /(I — z) = 1 for each z <E Z. 

The totally asymmetric exponential bricklayers process (TAEBLP) is ob- 
tained by taking 

(2.3) /(z) = e^(^-i/2). 

The construction of the bricklayers process with any nondecreasing / that is 
bounded by an exponential function is given in [B] on a set of tempered config- 
urations rj. This set consists of configurations with bounded asymptotic slope, 
the precise definition is given in 6^ . As certain desired semigroup properties are 
not fully proved, we avoid technical difficulties in the proofs of [8] by reproducing 
its results for the TAEBLP in the Appendix. 

2.2 The basic coupling 

We use a particularly simple form of the basic coupling which is made possible 
by the bricklayer representation: it is enough to define the structure of moves as 
described in jS] for a given side (left or right) of an individual bricklayer. Here 
is how to do it for a given bricklayer at site i. Given the present configurations 

uj} , oj^, . . . , oj" G fl, let m i-^ £(rn) be a permutation that orders the Ui values: 

^i < ^i , 1 < m < n. 

For simplicity, set 

p[m) : = /(t^f "•') and q{m) : = /(-wf "^), 

and the dummy variables p{0) = q{n + 1) =0. Recall that the function / 
is nondecreasing. Now the rule is that independently for each m — 1, . . . ,ri, 
at rate p{m) — p{m — 1), precisely bricklayers of a;^(™), u/i"^+^) ^ ...^ aj^(") 
place a brick on their right, and bricklayers of o;^^^-', w^'^-*, . . . , a;^(™~i) do not. 
Independently, at rate q{m) — q{m+ 1), precisely bricklayers oi u/^^\ w^*-^', . . . , 
u/(^n) pig^(3g £^ brick on their left, and bricklayers oi u/^™'^^\ u/(^n+2) ^ _ _ ^ u/(^) 
do not. Given the configurations to}, u^, . . . , u" £ ft, bricklayers at different 
sites perform the above steps independently. 

The combined effect of these joint rates creates the correct marginal rates, 
that is, the bricklayer of u;^(™) executes the move (j2.2[) with rate p{m) = 
f{u]^), and the same move on column /ii_i with rate q{'m) = /(— oj^ ). 

Notice also that, due to monotonicity of /, a jump of uf- without oJ' on the 
column [i, i+1] by the bricklayers at site i can only occur if /(w,f) < f{ujf) which 
implies wf > i^i- Also, a jump of w° without a;'' on the column [i — 1, i] by the 
bricklayers at site i can only occur if f{—Lu^) < /{—ujf) which implies ujf < ijj\. 
The result of any of these steps then cannot increase the number of discrepancies 
between the two processes, hence the name attractivity for monotonicity of /. 
In particular, a sitewise ordering ujf < uj\ Vi G Z is preserved by the basic 
coupling. 



The differences between two processes are called second class particles. Their 
number is nonincreasing. In particular, if, for processes w° and uj^ we have 
wf > cli^ for each i £ Z, then the second class particles between them are 
conserved. A special case that is of key importance to us is the situation where 
only one second class particle is present between two processes. 

2.3 Hydrodynamics and some exact identities 

From now on, we restrict our attention to the TAEBLP. Recall the jump rates 
(|2.2p . As described in [5;, the process has product translation-invariant station- 
ary distribution with marginals ^^ 

(2.4) /(z) 



/(z)! z{e) 

that turn out to be of discrete Gaussian type, see [1] for the explicit formula. 
P , E^, Var , Gov will refer to laws of a process evolving in this stationary 
distribution. The density g{9) : = E^(a;) S K is a strictly increasing function of 
the parameter S K, and can take on any real value by the Appendix of [5]. 
ji^, P^, E^, Var^, Cov^ will refer to laws of a density g stationary process. 
The hydrodynamic flux is 

n{Q) = W[f{u) + f{-uj)] = e''(^) + e-^(e). 

As / (|2.3p is convex and nonlinear, the Appendix of ^ applies and yields a 
convex hydrodynamic flux with 

(2.5) n"{Q) > 0. 

(This convexity property is quite natural, as H is the expected value of the 
convex rate function /. It has a nontrivial proof 5^ that uses total positivity 
or, a proof [7\ using correlation inequalities.) 
As in [Sj, we introduce 



oo 



*> °' z=y+l 

The Appendix of f5] applies to show that both ^^ and 'Jl^ are stochastically 
monotone in g. Denote by E the expectation w.r.t. the evolution of a pair 
(^~(')i ^(O) started with initial data (recaU (|2.1|) ) 



(2.7) ^-(0)-^(0)-^o^((g)A*^) 

and evolving under the basic coupling. This pair will always have a single 
second class particle whose position is denoted by Q{t). In other words, a;~(i) = 



w(i) — 5nu) ■ We reprove Corollaries 2.4 and 2.5 of "B^ in the Appendix that state 
that for any i £ Z and i > 0, 

(2.8) Var^(/i,(i)) ^ Var^(tj) • E|Q(i) - i| 

and 

E(g(i)) = F'^ • t, 

where T^^ — T-L'io) is the characteristic speed. Note in particular that in p.Sp 
the variances are taken in a stationary process, while the expectation of Q{t) is 
taken in the coupling with initial distribution (|2.7p . 

2.4 Results 

We repeat the results of [5], valid now for the TAEBLP. 

Theorem 2.1. Fix any density g E R, let the TAEBLP processes (^^(t), uj_{t)) 
evolve in basic coupling with initial distribution p.7p and let Q{t) be the position 
of the second class particle between w" {t) and w (i) . Then there is a constant 
Ci = Ci{g) G (0, cxj) such that for all 1 < m < 3, 

< hm mf — J- < hm sup — j- < 



Ci t^oo t2m/3 - j^^^ t2m/3 ^ _ ^u' 

where E on the right hand-side refers to the expectation in a coupled pair started 
from almost- equlibrium (|2.7I) . 

SuperdifFusivity of the second class particle is best seen with the choice 
m = 2: the variance of its position is of order t^^^. Next some corollaries. 
Notation [XJ stands for the lower integer part of X. 

Corollary 2.2 (Current variance). There is a constant Ci = Ci{g) > 0, such 
that 

1 Var^(Vf^j(i)) Var^(V^tjW) ^ r^ 

cT < H-S^ ws ^ 'TJ^^ ^75 < ^- 

Corollary 2.3 (Weak Law of Large Numbers for the second class particle). In 
a density-g stationary process, 

t 

Corollary 2.4 (Dependence of current on the initial configuration). For any 
V G M. and a > 1/3 the following limit holds in the L^ sense for a density-g 
stationary process: 

,2 g^ jjj^ ^LWJ jt) - hivti-iv,t\ (0) - tCHjg) - qH'{q)) ^ Q 



Recall that 



^ iu,{0), iiV<V^, 

i=lVt\-lVeti+l 

0, iiV = V^, 



hlVt\-lVet\{(^) 



E 



only depends on a finite segment of the initial configuration. Limit (|2.9p shows 
that on the diffusive time scale i^/^ only fluctuations from the initial distribution 
are visible: these fluctuations are translated rigidly at the characteristic speed 

Ve. 

Corollary 2.5 (Central Limit Theorem for the current). For any V ^ M. in a 
density-g stationary process 

hm ^-''^''l-^i^'^^ = Var^(c.) ■\V^-V\=: D, 

t^oo t 

and the Central Limit Theorem also holds: the centered and normalized height 
hiyt\ (t)/\/t ■ D converges in distribution to a standard normal. 

For convex rate zero range and bricklayers processes Corollaries 12.31 and 12.51 
were proved by M. Balazs [2j. 

The way we obtain the Theorem and the corollaries is simply proving that 
the assumptions formulated in J5| hold. For the sake of completeness, we repeat 
these, translated to our convex case. 

2.5 Microscopic convexity 

We start with the definition of microscopic convexity. This is just translated 
from the microscopic concavity property of 5 , where more detailed explanations 
and comments can be found. Convexity (j2.5p implies that the characteristic 
speed V^ = T-L'{q) is a nondecreasing function of the density g: 

(2.10) \< g =^ y^ < ye. 

The microscopic counterpart of a characteristic is the motion of a second class 
particle. Our key assumption that we term microscopic convexity is that the 
ordering (|2.10l) can also be realized at the particle level as an ordering between 
two second class particles introduced into two processes at densities A and g. 

Let A < £1 be two densities. Define /I^ + 1 as the measure that gives weight 
/ie(z — 1) to an integer z G Z (recall ()2.6p ). By the stochastic domination 
fl < /i^, we can let /i '^ be a coupling measure with marginals /I and /I^ + 1 
and with the property 

n^'^iy, z) : y < z} = L 



Let also fi^'^ be a coupling measure of site-marginals fi^ and /i^ of the invariant 
distributions, with 

(2.11) ^i^'S{{y,z)■.y<z} = l. 

Note the distinction that under fi^'^ the second coordinate is strictly above the 
first. 

To have notation for inhomogeneous product measures on Z^, let A — {Xi)i£z 
and g = {gi)iez denote sequences of density values, with A,; and gi assigned to 
site i. The product distribution with marginals 'j2^'>,so a,t the origin and ^^^^Bi 
at other sites is denoted by 

Measure Jl-'- gives probability one to the event 

{(7?(0), w(0)) : 770(0) < ojo{0), and 77,(0) < w,(0) for O^ieZ}. 

The initial configuration (77(0), io{0)) will always be assumed a member of this 

set, and the pair process (77(i), uit)) evolves in basic coupling. Notice that p,-'- 
is in general not stationary for this joint evolution. 

The discrepancies between these two processes are called the to — rj (second 
class) particles. The number of such particles at site i at time t is U!i(t) — rjiit). 
In the basic coupling the w — 77 particles are conserved, in the sense that none 
are created or annihilated. We label the uj — rj particles with integers, and let 
Xm(t) denote the position of particle m at time t. The initial labeling is chosen 
to satisfy 

■ • • < ^-i(O) < Xo(0) = < Xi(0) < • • • . 

We can specify that Xo(0) = because under /2-- there is an a; — 7; particle at 
site with probability 1. During the evolution we keep the positions Xi{t) of 
the a; — 77 particles ordered. To achieve this we stipulate that 

whenever an oj — 77 particle jumps from a site, 

(2.12) if the jump is to the right the highest label moves, 
and if the jump is to the left the lowest label moves. 

Here is the precise form of microscopic convexity for this paper. The assump- 
tion states that a certain joint construction of processes (that is, a coupling) 
can be performed for a range of densities in a neighborhood of a fixed density 
g. Recall (|2.ip for the definition of the configuration 6_. 



Assumption 2.6. Given a density p G K, there exists 70 > such that the 
following holds. For any A and g such that g ~ "fo < Xi < gi < g + ^q for 
all i e Z, a joint process (jjit), io(t), y{t), z{t))t>o can be constructed with the 
following properties. 



• Initially (77(0), io(0)) is Ji-'--distributed and the joint process {ri(-), a;(-)) 
evolves in basic coupling. 

• Processes y{-) and z{-) are integer-valued. Initially y(0) — z(0) — 0. With 
probability one 

(2.13) y{t) > z{t) for all t > 0. 

• Define the processes 

(2.14) ^"(t):-^(t)-5^^,„(t) and vt {t) ■ = ui^) + ^x^^,,(t)- 
Then both pairs (77, 77+) and {uj_^ , uj) evolve marginally in basic coupling. 

• For each 7 g (0, 70) and large enough i > there exists a probability 
distribution v^''^(t) on Z+ satisfying the tail bound 

(2.15) i^^'-'my ■.y>yo}< Ct^-W-^Vo'' 

for some fixed constants 3/2 < k < 3 and C < cx), and such that if 
g — j<Xi<gi<g + j for all i e Z, then we have the stochastic bounds 

(2.16) y{t) > -ve^''{t) and z{t) < vB^'^{t). 

Let us clarify some of the details in this assumption. 

Equation (|2.14p says that Q^{t) :— Xz[t){t) is the single second class par- 
ticle between 77 and 7;+, while Q{t) : — Xy(^t){t) is the one between cu" and uj_. 
The first three bullets say that it is possible to construct jointly four processes 
(7/, 77+, w^, Lo) with the specified initial conditions and so that each pair (77, tj), 
(77, 77+) and (w^, io) has the desired marginal distribution, and most importantly 
so that 

This is a consequence of (|2.13p because the uj — rj particles Xi (t) stay ordered. 

The tail bound (|2.15l) is formulated in this somewhat complicated fashion 
because this appears to be the weakest form our present proof allows. For the 
TAEBLP i'^''^ (t) will actually be a fixed geometric distribution. 

The assumptions made imply r]{t) < ui{t) a.s., and by (|2.14p 

77(i) < 77+ (i) < u{t) and 77(t) < ur{t) < ui{t) a.s. 

In our actual construction for the TAEBLP it turns out that while the triples 
(77, 77+, Lu) and (77, w^, cu) evolve also in basic coupling, the full joint evolution 
(77, 7;+, w^, Lo) does not. 

As already explained, the microscopic convexity idea is contained in in- 
equality (|2.13p . There is also a sense in which the tail bounds (|2.16p relate to 
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convexity of the flux. Consider the situation \i = X < g = Qi. We would expect 
the Lo ~ rj particle Xq{-) to have average and long-term velocity 

H{q)-H{\) 

R{X, Q) = , 

g-X 

the Rankine-Hugoniot or shock speed. By convexity T-L'{g) ~ V^ > R{X. g) > 
V^ = ■H'(A). A strict microscopic counterpart would be y{t) > > z{t). But 
this condition is overly restrictive. The distributional bounds (J2.16I) are the 
natural relaxations we use. 

Section [3] contains the proof of Assumption 12.61 for the TAEBLP. The proof 
of (|2.16p makes use of the particular exponential form (|2.3p of the rates. Un- 
fortunately, we do not have an argument for more general convex rates at the 
moment. 

There is one more assumption in [5] needed to state the main result. Con- 
stants C. , a. will not depend on time, but might depend on the density param- 
eter g, and their values can change from line to line. 

Assumption 2.7. Let (w^, uj) be a pair of processes in basic coupling, started 
from distribution (J2.7I) . with second class particle Q{t). Then there exist con- 
stants < aoj C < oo such that 

V{\Q{t)\>K}<C-^ 

whenever K > aot and t is large enough. 

Such an assumption is natural and easy to prove if the jump rates have 
bounded increments. Since / p.3p does not, this statement for the TAEBLP is 
nontrivial. We prove it in Section |4] for the TAEBLP. 

3 Proof of microscopic convexity 

In this section we verify that Assumption 12.61 can be satisfied. The task is to 
construct the processes y{t) and z(t) with the requisite properties. First let the 
processes (?](•), ^(•)) evolve in the basic coupling so that rii{t) < uJi{t) for all 
j G Z and t > 0. We consider as a background process this pair with the labeled 
and ordered w — 77 second class particles • • • < X_2(i) < X^i{t) < Xo{t) < 

Xi{t)<X2it)<--: 

At each time t > this background induces a partition {Aii{t)} of the label 
space Z into intervals indexed by sites i £ 1^, with partition intervals given by 

M^it)■.= {m : Xrnit)^l}. 

(For simplicity we assumed infinitely many second class particles in both direc- 
tions, but no problem arises in case we only have finitely many of them.) A4i{t) 
contains the labels of the second class particles that reside at site i at time t, and 
can be empty. The labels of the second class particles that are at the same site 
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as the one labeled m form the set Mx^it){t) = : {a™(t), a'"(t) + 1, . . . , 6"(t)}. 
The processes a"^ (i) and 6™ (i) are always well-defined and satisfy a™ (t) < m < 
b"'{t). Notice that 

(3.1) \Mx^it){t)\ = 6™(t) - a™(i) + 1 = ujx^it)it) ~ Vx^it){t). 

Let us clarify these notions by discussing the ways in which a™(t) and 6™(t) 
can change. 

• A second class particle jumps from site Xm{t~) — 1 to site Xm{t—). Then 
this one necessarily has label a'"(i— ) — 1, and it becomes the lowest labeled 
one at site X,„(i— ) = Xm{t) after the jump. Hence a™{t) — a™{t—) — 1. 

• A second class particle jumps from site Xm{t—) + 1 to site Xm{t—). Then 
this one necessarily has label b"^{t—) + l, and it becomes the highest labeled 
one at site X„i(t— ) = Xm{t) after the jump. Hence 6™(t) — b™{t—) + 1. 

• A second class particle, different from X„i, jumps from site A"„i(i— ) to 
site Xm{t—) + 1. Then this one is necessarily labeled 6™(i— ), and it leaves 
the site Xm(t-), hence 6"(i) = 6"(t-) - 1. 

• A second class particle, different from Xm, jumps from site Xm{t—) to site 
X,n{t—) — 1. Then this one is necessarily labeled a™(i— ), and it leaves 

the site Xm{t~), hence a™(t) = a"(f-) + 1. 

• The second class particle X„i is the highest labeled on its site, that is, 
m = 6'"'(t— ), and it jumps to site Xm{t—) + 1. Then this particle becomes 
the lowest labeled in the set Aix^(t-)+i = ■Mx„^{t), hence a"^{t) = m. In 
this case 6™(t) can be computed from (|3.ip . the number of second class 
particles at the site of Xm after the jump. 

• The second class particle X^. is the lowest labeled on its site, that is, 
TO = a™(t— ), and it jumps to site Xm{t—) — 1. Then this particle becomes 
the highest labeled in the set Mx,„(t-)-i = A^x„(t), hence 6™(i) — m. 
In this case a"* (t) can be computed from (|3.ip , the number of second class 
particles at the site of Xm after the jump. 

We fix initially j/(0) = z[Q) = 0. The evolution of (y, z) is superimposed 
on the background evolution (77, w, {Xm}) following the general rule below: 
Immediately after every move of the background process that involves the site 
where y resides before this move, y picks a new value from the labels on the 
site where it resides after the move. Thus y itself jumps only within partition 
intervals AAi. But y joins a new partition interval whenever it is the highest 
X-label on its site and its "carrier" particle Xy is forced to move to the next 
site on the right, or it is the lowest X-label on its site and its "carrier" particle 
Xy is forced to move to the next site on the left. 

These are the situations when y{t—) — 6'^(*^'(t— ) and at time t an ui — r] 
move from this site to the right happens, or y{t—) = a^^*^'{t—) and at time 
i an oj — 77 move from this site to the left happens. (Recall that the choice of 
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X-particle to move is determined by rule (|2.12p .') All this works for z in exactly 
the same way. 

Next we specify the probabilities that y and z use to refresh their values. 
Recall (j2.3p . To simplify notation, we abbreviate, given integers ry < w, 

(3.2) 

_ f{uj) - f{Lo - 1) _ fi-f]) - f{-T] - 1) _ e/^^"-") - e^('^-''-i) 



F\'h ^J - 


/(^) - /W 


f{-r^)-f(-u) 


qP(u-Vi) _ I 


and 








(3.3) 

q(77, u) = 


f{-u + l)-f{-u) 
f{-r,)-f{-w) 


fiv + l)-firi) 

/(^) - fiv) 


e^-1 


C/3(lj-7,) _ I ■ 



Notice that both p{ri, uj) and q{r], uj) only depend on uj — rj. Therefore, with a 
little abuse of notation, we write p{io — ri): = p(ri, lu), q{uj — rj) : = q(ri, lu). Then 

p(l) = <7(1) = 1, p{d)>q{d), p{d)+q{d)<l for 2 < d G Z. 

When y and z reside at separate sites, they refresh independently. When they 
are together in the same partition interval, they use the joint distribution in the 
third bullet below. 

• Whenever any change occurs in either a; or 77 at site Xy(^t-){t~) and, as a 
result of the jump, a^(*^^(i) 7^ a^'^*^^(t), that is, y{t~) and z(t~) belong 
to different parts after the jump, we abbreviate 

P = PiVx.^t-) it) (*)> '^x,^,_) it) (t)) q = 9fe,(._, (t) (i), '^x,^t-) («) (0) 



of (|3.2p and p.3p in the formulas below. These depend on the values of 
the respective processes at the site where the label y can be found right 
after the jump. In this case, independently of everything else. 



(3.4) yit) 



a^(*")(i), with prob. g, 

by'-*-\t) - 1, with prob. l-p-q, 
by^*^\t), with prob. p, 

except for y[t) : — ay^^^\t) — 6^(*~^(t) when the difference 

'^x„(,_,(t)(i)-»7x„(,_,(t)(i) 

is 1. Notice that the second line in (13.41) has probability zero when this 
difference is 2. 

Whenever any change occurs in either w or 77 at site X ^(i_\{t—) and, as a 
result of the jump, a^*^*^^(i) ^ a^'*^^(t), that is, 2/(t~) and z{t—) belong 
to different parts after the jump, we abbreviate 

P = p('7x,(,_)(t)(i), a;x,(,_)(t)(i)) q ^ q{f^x,^,_^(t){t), ^x,^,_^(t){t)) 
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of p.2p and p.3p in the formulas below. These depend on the values of 
the respective processes at the site where the label z can be found right 
after the jump. In this case, independently of everything else, 



(3.5) 



z{t) 



z(t-) 
z(t-) 



except for z{€) 



(t), with prob. p, 

(t) + 1, with prob. 1 — p — g, 
6"^(*~)(i), with prob. 9, 

a^(*-)(f) — 6^(*^)(i) when the difference 



•^x. 



,(t)W-'7Jf,(,_)(t)W 



is 1. Notice that the second line in p.Sp has probability zero when this 
difference is 2. 

• Whenever any change occurs in either w or 77 at sites Xy{t-)^~) or 
^z(t~){^—) and, asaresult ofthe jump, a^^*~-'(t) = a^'*^)(t), that is, y(t—) 
and z(t—) belong to the same part after the jump, that is, XyU\(i) = 
X^u\(t) then we have 



and we abbreviate 

P = PiVXy^,_)(t) (0, '^X,;,^) (t) {t)) 



and ??x„(,_)(t)(i) = ??x,(,_)(t)(t), 



q{Vx,f,_^{t){t),^Xy^,_^{t){t)) 



of (|3.2I) and (|3.3p in the formulas below. These depend on the values of 
the respective processes at the site where both the labels y and z can be 
found right after the jump. In this case, independently of everything else. 



(3.6) 



z{t) 






, with prob. q, 



r,v(t-)u\ ) ' ^i*h P'"°b- (P - 9) A (1 - p - q), 



Kt) 



: = < 



6!^(*-)(t) 
it) 



-,y{t-) 



by^*-\t) 



, with prob. [2p — 1] 
1 



,a(t-) 



, with prob. [1 — 2p]^, 

, with prob. {p — q) A {I — p — q), 



ty^'-Ht) 
by^'-^t) 



, with prob. q. 



except for y{t) = z{t) : = a«(*-)(^) = by^'^-\t) when the difference 

^ Xy,t^)(t){i) ~ 'nXyn-)(t)^^) i^ 1- Notice that the second, the fourth and the 

fifth lines have probability zero when this difference is 2. 
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The above moves for y and z always occur within labels at a given site. This 
determines whether the particle Q{t) : = Xy(^t){t) or Q''(i) : = Xz{^t){t) is the one 
to jump if the next move out of the site is an w — rj move. 

We prove that the above construction has the properties required in As- 
sumption 12.61 First note that the refreshing rule p.6|) marginally gives the 
same moves and probabilities as p.4p or (|3.5p for y{-) or z(-), respectively. 

Lemma 3.1. The pair (w^, uj) :— {ui~ 5_x , w) obeys basic coupling, as does 
the pair {ri, ^+) : = (?7, »7 + ^xj- 

Proof. We write the proof for (w^, uj). We need to show that, given the con- 
figuration (?7, w, {Xm}, y), the jump rates of (w~, w) are the ones prescribed 
in basic coupling (Section 12. 2|) and by p.2|) . As mentioned in Section [2?2| the 
effect of bricklayers determine the evolution of processes. Notice first that an 
u ~ T] particle can only jump away from a site « if a bricklayer of w or 77 moves. 
As the moves (|3.4I) or (|3.6p by themselves never result in a change of X^ (.)(•), 
any move of Q from a site i is a result of a bricklayer's move at site i. Therefore, 
we see that moves initiated by bricklayers of w at sites i ^ Q happen as well to 
w", as required by the basic coupling. The only point to consider is moves by 
the bricklayers at site i = Q. We start with them putting a brick on their right. 
Since the last time any change occurred at site i, y chose values according to 
(|3.4I) or (13. 6p . Notice that (|3.4p and (13. 6p give the same marginal probabilities 
for this choice. Hence 



(3.7) 

y took a value less than b^ with probability 1 — p — 

and 

(3.8) y took on value b^ with probability p - 



f{u, - 1) - fjrji) 



!{lO,) - fJLU, - 1) 



as given in p.4p . Notice that p.7p happens with probability zero if oji = r/i + 
1. According to the basic coupling of rj and w, the following right moves of 
bricklayers at i can occur: 

• With rate /(wi) — f{rii), w jumps without 77. The highest labeled second 
class particle, Xbv jumps from site i to site i + 1. 



With probability (|3.8p Xy = Q jumps with X^y . In this case 

ujr{t-) = uj.,{t-) - 1 = uj.,{t) - ujr{t) 
since the difference Q disappears from site i. Also, 

w,^i(i-) = w,+i(i-) = UJ^+l{t) - 1 = W,^i(t), 
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since the difference Q appears at site i+l. So in this case cj undergoes 
a jump but uj_~ does not, and the rate is 

[/(-.) - /(..)] • y(^^) _ ;(,^) - /(-^) - /K ). 

— With probability (|3.7p Xy = Q does not jump witlr Xi,y , since it has 
label less than b^ (this probability is zero \i ijji — rji + l). In this case 
ui~ and oj perform the same jump and it occurs with rate 

[/(-o - Km)] ■ ^^^; ~!^ ~J^f - /(-D - fim)- 

JK^i) -Jim) 

• With rate /{rji), both bricklayers of r/ and u at site i move. No change 
occurs in the w — 77 particles, hence no change occurs in Q. This implies 
that the process w~ jumps as well. 

Summarizing we see that the rate for the bricklayers of (w~ , w) at site i to lay 
brick on their rights together is f{u~), and the rate for the one of w to move 
without w^ is /(wi) — /(aj~). This is exactly what basic coupling requires. 

Consider now bricklayers at site i = Q putting a brick on their left. Since 
the last time any change occurred at site i, y chose values according to p.4p or 
(13:61). Hence 



(3.9) y took on value a^ with probability q — — -j — 

I{-ni) - I{-^i) 

and 

(3.10) 

y took a value higher than a^ with probability 1 — q = — — r — 7 — , 

f{-m) - f{-'^i) 

as given in (|3.4I) . Notice that (I3.10p happens with probability zero if w^ = ?7j + 1. 
According to the basic coupling of r/ and w, the following left moves of bricklayers 
at i can occur: 

• With rate f{—rii) — f{—LOi), rj jumps without to. The lowest labeled second 
class particle, Xav jumps from site i to site i — 1. 

— With probability (|3.9p Xy — Q jumps with Xav ■ In this case 

u;~{t~) = uj,{t-) - 1 = Lu,{t) - 1 = Lor{t) - 1 

since the difference Q disappears from site i. Also, 

^i-i{^-) = ^i-iii-) = ^i-i(i) = ^i-iii) + I7 

since the difference Q appears at site i + l. So in this case w" 
undergoes a jump but w docs not, and the rate is 
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— With probability p.lOp Xy = Q does not jump with Xav, since it 
has label more than a^ (this probability is zero if Wi = 77^ + 1). In 
this case none of lu~ or to move; this occurs with rate 

Uf \ ft M /(-^O - /(-'^l + 1) ff ^ f, -X 

m-n.) - /(-.)] • ;(_,^) _ ;(^^) = /(-'^) - /(-. )• 

• With rate /(— Wi), both bricklayers of rj and a; at site i move. No change 
occurs in the uj — i] particles, hence no change occurs in Q. This implies 
that the process w~ jumps as well. 

Summarizing we see that the rate for the bricklayers of (w~ , w) at site i to lay 
brick on their rights together is f{—uji), and the rate for the one of w" to move 
without w is f{—uj^) — f(~LjJi). This is exactly what basic coupling requires. 
A very similar argument can be repeated for (rj, r]~^). D 

Lemma 3.2. Inequality (j2.13p y > z holds in the above construction. 

Proof. Since no jump of y or z moves one of them into a new partition interval, 
the only situation that can jeopardize p.l3|) is the simultaneous refreshing of y 
and z in a common partition interval. But this case is governed by step p.6p 
which by definition ensures that y > z. (When 6^^*~''(i) = a*'(*~)(i) + 1, we 
have, by dSH]), Wx„(t_)(t)(0-'7x„(t_,(t)(0 = 2, and hence p of dMI) is more than 
1/2. Therefore the probability of the step in line 4 of (13.61) is zero.) D 

Define the geometric distribution 

-e"''), m>0 




(3.11) ^^(m) , ^ 

m < 0. 

d 

Lemma 3.3. Conditioned on the process {rj, lo), the bounds y(t) > —v and 

d ~ 

z{t) < v hold for all t > 0. 

To avoid unnecessary complications with negative values, we show the proof 
for z{t). Notice that both the statement and the behavior of y{t) is reflected 
compared to z{t), hence the proof is the same for the two processes. The 
argument consists of three steps. 

Lemma 3.4. Let Z be a random variable with distribution u, and fix integers 
a < b and rj < uj so that uj — rj — b — a + 1. Apply the following operation to Z : 

(i) if a < Z < b, apply the probabilities from p.Sp with parameters a, b, rj, w 
to pick a new value for Z ; 

(ii) if Z < a or Z > b then do not change Z . 

Then the resulting distribution v* is stochastically dominated by v . 
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Proof. There is nothing to prove when b — a, hence we assume b > a or, 
equivalently, w — 77 = 6 — a + l>2. It is also clear that f* (m) = v{m) for m < a 
OT m > b. We need to prove, in view of the distribution functions, 

m m 

for all a < TO < 6. Notice that v* gives zero weight on values a + 1 < rra < 6 (if 
any), and also that the display becomes an equality if m = 6. Therefore, it is 
enough to prove the inequality for in = a: 

(3.12) y*{a)>v{a), 
and m = b — 1: 

h-l fc-l 

(3.13) ^v*{^)>^v{i) that is, v*(b)<v{b). 
e=a e=a 

Notice that (|3.12p is trivially true for a < 0. For a > we start with rewriting 
the left hand-side of (|3.12p with the use of p.5p . (|3.2I) . and the abbreviation 
d = cj — ry = 6 — a+1: 



\a)^pid)-j2m 



Q0d _ g/3(d-l) 



g,3d _ I 



fSa _ p-/J(6+l) 



(e-'^'^ - e 



= --'"■ J_, -{1-0-^') = Ha). 
As for p.l3p . both sides become zero if 6 < 0. For 6 > we have 

6 
l.*{b)^q{d}-Y,,.{i) 
t=a 

< ^^"-^ . (e-f^" - e-^(''+i)) 
- el3d _ I ^^ ^ I 



= e-P' 



e 



3/3d _ I 



(e/5('^-i)_e-^) = K&). 



Lemma 3.5. The dynamics defined by ()3.5p is attractive. 



U 



Proof. Following the same realizations of p.Sp . we see that two copies of z(-) 
under a common environment can be coupled so that whenever they get to the 
same part Aii, they move together from that moment. D 
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Proof of Lemma \3.S[ Initially z(0) = by definition, which is clearly a distri- 
bution dominated by i^ of (|3.11l) . Now we argue recursively: by time t the 
distribution of z{t) was a.s. only influenced by finitely many jumps of the en- 
vironment, which resulted in distributions vi, then 1/2, then v^, etc. Suppose 

d 

Uk < ly, and let v* be the distribution that would result from v by the k + V^ 

d d d 

jump. Then i^k+i < i^* by Vk < i^ and Lemma 13.51 while ly* < v hy Lemma 
[3:1 D 



4 A tail bound for the second class particle 

In this section we prove that Assumption 12.71 holds for the TAEBLP model. 
The difficulty comes from the fact that jump rates of the second class particle, 
being the increments of the growth rates (|2.3I) , are unbounded. First recall the 
coupling measure /i^'^ of (|2.1ip and notice that it gives weight one on pairs of 
the form {y, y) if A = g. Define also ^^^'^°'^^e by 

I 0, otherwise. 

With these marginals we define the shock product distribution 
(4.1) if^"""^ '^ — Ai''+^'''+^ • m'''""'' '^ • m'''^ 

i<0 i=0 i>0 

a measure on a pair of coupled processes with a single second class particle at 
the origin. 

Lemma 4.1. The first marginal oj^s'i-ocfce jg ^^g pj-Qf^^Qi distribution 



i<0 i>0 



while the second marginal is 



(4.2) (g)/i^+^-(g)M'- 

Proof. The first part of the statement and the second part, apart from i — Q, 
follow from the definitions. The nontrivial part is 

valid for the second marginal at i = 0. This is specific to the definition (12. 4p of 
/u^, and of the exponential rates (|2.3p . and to prove it wc write, with 9 = 0{g), 

j(^) e«- 1 e(«+^)- 1 

^e(^z - 1) = — 



/(z)! Z{0) f(z)l e«+'3/2Z(6l) ■ 
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Summing this up for all z S Z gives one on the left hand-side, hence leads to 

2(e+/3) 






which also implies 



q{o + /?) = ^ iog(^(^ + P)) = 9{e) + 1. 



We conclude that 



p(e+/3)z 1 

which finishes the proof of the lemma. □ 

The translation of ^^^hock g jg (denoted by 

i<ik i—k i>k 

The main tool we use is Theorem 1 from [3], which we reformulate here. fJ,S{t) 
will just denote the time evolution of a measure fi under the process dynamics: 

Theorem 4.2. In the sense of bounded test functions onfl x fl, 

+ (e"^(e) _ e-0(e+i)) . (rfe_i^"''°^''-e - /j,«''°^*=e)_ 

The first interesting consequence of this theorem is that the measure ^^^^""^^ e 
on a coupled pair evolves into a linear combination of its shifted versions. Sec- 
ond, notice that (|4.3p is the Kolmogorov equation for an asymmetric simple 
random walk. Indeed, this theorem implies the following 

Corollary 4.3. Let the pair (^^(0), ^(0)) have initial distribution ^'^^°<^^e gjg_ 
fined by (|4.1|) . Then its later distribution evolves into a linear combination of 
translated versions of fj,"^"'^'' ^ : at time t the pair {£,^{t), £,(0)) has distribution 

oo 
fc— — oo 

where Pk{t) is the transition probability at time t from the origin to k of a 
continuous time asymmetric simple random walk with jump rates 

e«(e+i) _ e«(e) to the right and e~^^e) _ Q-~e(e+i) ^^ ^^g i^j^^ 

In particular, Q^{), started from an environment ^'^^"'^^ ^ , is a continuous time 
asymmetric simple random walk with these rates. 



20 



Although the corollary is quite natural, let us give a formal proof here. First 
some notation. {^~{-), ^(•)) will denote a pair of processes evolving under the 
basic coupling, g will be a bounded function on the path space of such a pair, 
and for shortness we introduce Qt for the whole random path, shifted to time t: 
&t = (C~ (*+•)> i(*+'))- Expectation of the process, started from TkfJ^^"'^^^, will 
be denoted by E^*^'. Notice that under E^*^) we a.s. have a single position Q^{t) 
where the coupled pair differ by one, this is the position of the single conserved 
second class particle. With some abuse of notation we also use E^i '^' for the 
evolution of the pair {^~{-), ^i-)), started from the specific initial state (^~, ^). 

We aim for proving the semigroup property of S{-). The first step is 

Lemma 4.4. Given times < s < t and fc G Z, 

E(")[5(et)|g«(s) = fc] = EW[5(e,_,)]. 

Proof. The left hand-side is 

p(o){Qe(s) = k} 

_ EW[E(r(^).i(-))g(9t-.);Q^(.)^fc] 
p(o){Qe(s) = k} 

E p(°HQ«W -.?}EO)[Er(o).i(o))g(e,_,) ; g«(o) = k] 

^ iez 

P(o){Q«(s) = k} 

_ p(o){g«(s) = fc}EW[E(r(o).i(o))g(et_,) ; g«(o) = k] 

" P(0){QC(s) = fc} 

= EW[.9(e*-.)], 

where in the second equality we used that the distribution at time s is a linear 
combination of shifted versions of pt^^hocke q 

Next we prove the Markov property for Q^{-). 

Lemma 4.5. Let n > be an integer, (pi, i — 0, . . . , n bounded functions on 
Z, and ~ to < ti < ■ ■ ■ < tn- Then 

n n 

^^''^Xl^^{QHti)-QHU-l))=Y{'E.^''^V^{Q^{t^-t,-^)). 
i=l i=l 

Proof. The statement is trivially true for n = 1. We proceed by induction, and 
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assume the statement is true for n — 1. Then 

n 

n 

= ^p(o){g«(ii) = j}<^i(j) • E(o) [H^.iQHu) - qHu-i)) I QHti) = j 

n 
n 
n 
n 

=nE(°v.(Q^(i.-t.-i))- 

i=l 

The second equality uses Lemma [4.41 the third one uses the fact that 0's only 
depend on Q^-differences, and the fourth one follows from the induction hy- 
pothesis, n 



Proof of Corollary \4-3\ We know that at any fixed time i > the distribution 
of {£,~ {t), ^{t)) is a linear combination of shifted versions of ^^hockg^ The shift 
is traced by the second class particle Q^{t), therefore the differential equation 

(4.4) ^p(0){Q«(t) = fc} 

= {Ss+D _ c^Ce)) . (p(o){g«(t) = fc + 1} - p("){Q«(i) = k}) 

+ (c-«(^°) - c-'^s+i)^ ■ (pW{Q«(i) = fc - 1} - p(o){g«(t) = fc}) 

follows from (|4.3p . In the above lemmas, we also proved that Q^{t) is Marko- 
vian (annealed w.r.t. the initial distribution of {£,~ , f)). As there exists only 
one Markovian process with Kolmogorov equation (|4.4p of the simple asymmet- 
ric random walk, we conclude that the process Q^(-) with initial environment 
^shock e jg g^j^ asymmetric simple random walk with rates as stated in the Corol- 
lary, n 

Lemma 4.6. Let (w^, w) be a pair of processes in basic coupling, started from 
distribution ()2.7p . with second class particle Q{t). Then there exist constants 
< ao, C < oo such that 

P{\Q{t)\ > K} < e-^^ 
whenever K > aot and t is large enough. 
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Notice that this imphes that Assumption 12 . 71 holds for the TAEBLP. 

Proof. The proof uses auxihary processes to connect the above arguments to 
the setting of Assumption 12.71 Define the pair 



\ig,g+l), fori<0, 

(Ai, g,) := -> f ■ ^ n 

I [g, g), tor t > Q. 

Draw the pair (C(0), C(0)) from the product distribution of coupling measures 

(Em 



Then ^(0) has distribution 



IM 



i<0 i<0 



in agreement with (I4.2p . 

Let now the pair (C('): ?(')) evolve in the basic coupling, and let them play 
the role of (?/(■), w(-)) of Section[3l This results in the pair {({■), C'''(')) ^ith a 
second class particle Q''{-) and the pair {^~{-), £,{■)) with a second class particle 
Q^{-) such that Q''{t) < Q^{t), see Lemma 13.21 Therefore the random walk 
result in Corollarv 14.31 on Q^(-) yields the desired estimate for Q''{t). Finally, 
notice that the distribution of aj~(0) in Assumption 12 . 71 and of C(0) above only 
differ by 0Jq{0) ~ /2^, while Co(0) ~ M^- Therefore 

P{Q{t)>K}= ^ P{Q(t)>K\co-{0) = z}-f,^{z)i[i^y 

Z — — 00 r" \ ) 

Z— — 00 r' \ J 

< [ Y: P{Q^ W > K I Co(0) = z} ■ M^(z)] ^ [ ^ ^ ^^^ 



,=-oo ^''^y) 



p 



oo 



ti'iyf 



y=-oo f'(y) 



We are done as soon as we show that jl^ (y) / fJ-^ (y) is uniformly bounded in y. 
With the exponential rates (|2.3p one obtains from 

^ '-y-' z=y+l fc=l 
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This is uniformly bounded for large j/'s since then ye ^^ < \. For large negative 
y's one uses the equivalent form 



of (|2.6p and writes 

C J2 ig~z)e-^^^-7r+i(y-fr^cJ2ik-y + g)e-^ 



V'^^iv) _r^ V- /„ .^„-l(^-l)' + 4(w-4)' _^V-/,, .. , .^„-*fc^+/3fca-efc 



which is again uniformly bounded for large negative y values. 
To show a lower bound on Q{t)^ start with 

{(e, e), for z < 0, 

{q- 1, e), for i = 0, 
(e, e-1), for i > 0, 

and the coupled pair (C(0), C(0)) in distribution 

Now the roles of the pair (C(-)i C^(')) with a second class particle (3^(-) and 
the pair (C (Oi C(')) with a second class particle Q^{-) are interchanged and we 
have Q'^{t) > Q^{t). The random walk estimate on Q^ and a Radon- Nikodym 
estimate similar to the one above completes the proof of the lower bound. D 

Appendices 

A Covariance identities for bricklayer process 
with exponential rates 

The purpose of this appendix is to prove the variance formula for stationary BLP 
under the following exponential bound assumption on rates: for some C, /3 < cx), 

Assume the height process is normalized initially by ho{0) = 0. 

Theorem A.l. Fix z € Z. In the stationary infinite volume process with time- 
marginal distribution a;(i) ^ fi^ , 

(A.l) Var[/i,(i)] =^|n-z|CovK(i),Wo(0)]. 
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Formula (|A.1[) was already proved in [5] for a general class of processes 
that includes ZRP and BLP. However, the proofs in [8] were carried out under 
the assumption that certain semigroup and generator calculations work. This 
presents no problem when the single-site state space is compact (such as exclu- 
sion processes ) for then one has a strongly continuous semigroup on the Banach 
space of continuous functions on the state space of the process. For BLP with 
superlinear rates, only some rudimentary features of the usual semigroup pic- 
ture have been established in [6] . Hence the need to justify (lA.ll) . We use the 
finite- volume auxiliary processes as introduced in [B] . 

To prove Theorem lA.ll we show that the infinite- volume stationary process 
is a limit of finite- volume {£, r, 9) processes as —£, r — >■ oo. A preliminary form 
of (|A.ip is true for an (£, r, 9) process by simple counting. (See (|A.9p below 
and its expanded form on lines (jA.lip - (jA.14[) .') The technical work goes into 
establishing moment bounds that are uniform over £ < < v. These in turn 
permit us to take the —£, r — ?► oo limit in the proto- formula (jA.lll) - (jA.14p . 

This appendix is based on the construction of the infinite-volume BLP h{t) 
given in [B]. Article |^ utilized two types of finite- volume processes: the [£, r] 
processes whose height variables were denoted by h}^'^'{t), and the {£, r, 9) pro- 
cesses with height variables denoted by g^^'^'^\t). 

The [£, r] evolution is a straightforward restriction of the full system into the 
finite interval [£, r], with generator 

t-i 

This generator defines a countable state space Markov process that evolves over 
the sites £, . . . , r: the jump oj -^ uA^''^+'^) happens at rate /(wi) -I- /(— oji+i), 
independently at different sites i, but only for £ < i < r — 1. Columns outside 
the interval [£, r] are frozen for all time. The virtue of this process is monotone 
dependence on the interval [£, v\. In [6] the infinite- volume process h(-) was 
defined as the a.s. increasing limit of the processes h} ''^'{t) as — ^, r — > oo. 

The {£, r, 6) process has also the correct boundary currents that make the 
i.i.d. product measures /i^ invariant for the finite collection of increment vari- 
ables {Wj- ''^' . ^ < i < r}. The generator is 

r-l 

For a concrete construction of the processes, we imagine that bricklayer at 
site i has two unit rate Poisson processes N^ and N> on the first quadrant 
R^ of the plane. These govern his brick-laying action to the left (L) and right 

(i?). A Poisson point {t,y) in N> ' such that y < f{—uji{t—)) signals a brick 
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to be laid on [i — l,i] at time t, while a Poisson point {t,y) in iV| ' such that 
y < f{^i{t~)) signals a brick to be laid on [i,i + 1] at time t. Shift of these 
planar Poisson processes by time t will be denoted by StN. 

The construction of the finite systems in terms of these Poisson processes 
provides the usual jump chain-holding time construction of a continuous time 
Markov chain on a countable state space. After each jump, the holding time 
and the next state are read from the Poisson processes. By the strong Markov 
property of the Poisson processes this is equivalent to looking at freshly sampled 
exponential variables with appropriate rates. In both [£, r] processes and {£, r, 9) 
processes the increment H{t) — H{0) of the maximal height 

(A.2) i/(i)= max /i,(t) 

«— l<j<r 

is bounded by a Poisson process. Hence explosions do not happen [6, Sect. 3.1- 
3.2]. 

In the next section we show that the {i, r, 9) processes converge to the 
infinite volume stationary process. Then we develop moment bounds through 
martingales, uniformly in £ < < r. After Section lA.l I wc drop the superscripts 
[£, r] and {£, r, 9) to ease notation. 

A.l Convergence of (£, r, 6) processes 

The infinite-volume process h{-) is defined as the a.s. increasing limit of the 
processes h}^''^'{-) [SI Theorem 2.2]. Lemma 7.1 in [B] shows that if the initial 
state w(0) is /i^-distributed then the increment process lu{-) is stationary. 

Calculations in this appendix are done in a stationary finite- volume (£, r, 9) 
process. Hence we need to show that this process also converges as —i, r — ?> oo 
to the stationary infinite-volume process. 

It will be convenient to represent the processes by measurable mappings of 
the initial configuration h and the Poisson clocks N_ on (0, t]: 

^(^.^^»)(t) = ^(^^^'«)(/i,iV), h^^''\t)^<^f''\h,N) and h{t) ^ '^t{h,N). 

Then the construction of the process h{t) in P; can be expressed as follows: for 
any initial hCzQ, any m,T < oo, 

(A.S) $t,i(/i, N) = $5 '1 {h, N) for -m<i<ma.ndO<t<T 

for all large enough —£, r, almost surely. 

Lemma A.2. Let h(t) be the infinite-volume process whose increment process 
tu{t) is stationary with marginal ui{t) ^ pfi . As — ^,r — >■ oo, (^(^' '^'^)(.) — ^ /i(.) in 
the following sense: given any m,T < oo, 

gf '"■'^^i) = h^{t) for -m < i < m and < t <T 

for all large enough — ^, r. 
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Proof. The proof of Lemma 7.1 on p. 1243 of [B] shows that there exists a 
(nonrandom) time to = to{9) > such that, for any m, gl ''^' '(t) = /i[ ' (t) for 
—TO < i < m and < i < io if — ^ and r are large enough. Combined with (jA.Sp 
we have the statement on the time interval [0,to]. 

To get to an arbitrary time T, the claim is proved by induction up to time 
kto > T, for a positive integer k. Since process h{-) stays in ft, (|A.3[) can be 
applied to restarted processes. 

The other ingredient of the induction step is control of discrepancies between 
processes 5, (•) and h\ ' (•). For this purpose let B be the event that in 

the process {g^°'''''^''{kto + t) : < t < to} either all columns in the range 
{ [M/4\ , . . . , [Z/2J + 1} grew, or all columns in the range { [c/2j , . . . , [3r/4] } 
grew. By [6l Corollary 5.5], P(/B) < Cie~^^^^^^^'^\ and this bound is independent 
of a, b. 

Discrepancies originate at the edges £ and c. On the event B"^ the (a, b, 9) and 
[£, r] processes started from a common initial configuration are indistinguishable 
inside (£/2,r/2) because discrepancies have not had a chance to penetrate into 
{i/2, r/2). By the monotonicity properties of these processes the columns of the 
(£, r, 9) process never go below those of the [£, r] process [U Lemma 3.2]. Hence 
it is enough to use the event B"^ to suppress growth in the {i, r, 9) process. We 
leave the details of this induction argument to the reader. D 

A. 2 Martingales 

Since the rates are unbounded, we begin by stating a general lemma about 
countable Markov chains. Let S* be a countable state space, Q a generator 
matrix, q^ = —qx,x — J2yy^x l^.v ^^^ total rate to jump away from state x G S. 
Let P(t) be the semigroup associated to Q, in other words the minimal positive 
solution of the backward equation P'{t) = QP, -P(O) = / [5]. The next lemma 
is proved with standard techniques. 

Lemma A. 3. Let v he an initial distribution on S and T < cxd. Assume 

T T 

(A.4) j E^lqx^ds^ I ^v{x)^px,y{s)qy ds < oo. 

^ y 

Let ip be a bounded function on S. Then under P'^ , for t G [0,T], the process 
a{t) = L Q(p{Xg) ds is well-defined and in L^{P'^), and the following process is 
a martingale: 

t 

(A.5) Mt - ^(Xt) - ifiiXo) - [ Qv{X,) ds. 



Now apply this to the [^, t] and (^, r, 9) processes. To simplify notation, let 
f.{ui{s)) = f{u,{s)) + f{-u:,+,{s)) 
be the rate of growth of the column height hi{s). 
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Lemma A. 4. Fix £ < < x. Consider either the [£, r] process or the {£, r, 9) 

process and in either case denote the height variables by hi(t). Let v be an initial 
distribution such that for some c > /3 

(A.6) ^"(e'l'^ol + '^-^l'^'l) <oo. 

Then for any 1 < p < oo and any index i this process is a martingale: 



(A.7) M, = h^{t) - /.f (0) - / f.Ms)){{K{s) + ir - hUs)) ds. 



Proof. In the [£, r] process the total jump rate out of state h is 



r-l 



qh 



'/L=E/^(^)^2^^^'"''^2(r-£+l)e''^.-^ 



I". I 



In the {£, c, 6) process the total jump rate out of state g is 

Qt = % + fi-^i) + /K) + £" + e-o <2{x-£+ l)e''S'=^l"'l +6^ + 6" 
In the [£, r] process under a fixed initial configuration h, 



W 



sup e 

se[o,T] 



/5EL.I'^.(s)l 



<e^i:i=A-',\-E[e^y(T)^ 



where Y{-) is a Poisson process of rate A = 2/(0)(c — £). This comes from the 
observation that the process 

r r 

z;(i) = ^|...(i)|-^N(0)| 

increases only when a local maximum column grows, and this happens at rate 
at most 2/(0) (r — £). Then, under the initial distribution v, 



E'' 



(A. 



sup g^(s) 

sG[0,T] 



< w 



sup e 

sG[0,T] 



/J EL. I", (s) I 



< E''(e'^^'=^l'^'l) • exp{A(e'^^ - 1)} < oo. 



In the (£, t, 9) process the same idea works except the rate for the process 
that bounds v{t) is A^ = 2/(0)(r - ^ + 1) + e^ + e-**. 

Let F{x) — (6 A x) V (—6) be a truncation function. Now that assumption 
(jA.4[) has been verified, for any integer < 6 < oo and ^ < i < c — 1 (JA.SP gives 
the martingale 



mI'^ ^ F{h,{t)Y - F{hmY - f.{ui{s)){F{h,{s) + ir-F{h,{s)r)ds. 
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So for an event A ^ F^ and s < t 



As 6 /* oo dominated convergence takes each term to the desired hmit. This is 
justified by the fohowing. Restrict to an interval s,t € [0,T]. The rate in the 
last expectation is bounded as in 

h{u{u)) < sup e'5(l".WI+l-'.+iMI) 

selo,T] 

and the random variable on the right has a finite L^-norm for some p > 1 hy 
a bound of the type in (jA.SP and the assumption that c > /? in (jA.6| . For the 
height we have 

Mt)\ = h+it) + h-{t) < H{t) + V(0) < {H{t) - H{Oj) + H{0) + M0)\ 

where H{t) is the maximal height of (jA.2[) . The increment H{t) — H{0) is 
stochastically bounded by a Poisson process while H{0) and |/ii(0)| have all 
moments by assumption (jA.6|) . We conclude that (|A.7|) is a martingale. D 

A. 3 Bounds for the {i, r, 6) process 

Henceforth consider stationary {£, r, 0) processes with /i^ marginals and initial 
height normalized by /io(0) — 0. The increment process is denoted by Lu{t) and 
the height process by h(t). In this process the column heights hi{t) ioT i < £ — 2 
and i > X + 1 are frozen at their initial values. We start with moment bounds 
that hold uniformly in £ < < r. 

Lemma A. 5. Fix 6. For I < p < oo there is a constant C = C{p, 9) such that 
in all stationary {£, r, 9) processes with marginal distribution pP , for all t > 
and i G Z, 

E[ih,it)-h,iO))P] <e". 

The bound is valid also for the boundary columns i = £ — I and i = r, and also 
for the infinite-volume stationary process with marginal distribution jjfi . 

Proof. Abbreviate hi{t) = hi{t) — hi{0). It suffices to consider the {£, r, 9) 
processes because the bound extends to the —£, r — >■ oo limit by Fatou's lemma 
and Lemma I A. 2 1 The increments U!i have all exponential moments under /x so 
assumption (|A.6p is satisfied. In particular, hi(t) has all moments. 

Fix £ < < c. Columns hi for i ^ [£ — l,r] are frozen at their initial values 
and need no argument. Consider the case £ < i < t — 1. It suffices to consider 
a positive integer p. In the next calculation, use martingale (jA.7|) , the fact that 
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ft-i(s) is a nonnegative integer to bound a lower power with a higher one, and 
apply Holder's inequality: 



* p-i 



/c=0 ^ ^ 

t 
<CpjE[Moi{s)){lvh{s)Y-']ds 



t 



t 
< Cp,g f {-Ehisf + I) ds. 



Rewrite this as 

t 

Eh,{ty + i<i + Cp,8 l{Ehi{3) + 1) ds 



and now Gronwall's inequality gives the conclusion. 

The boundary columns i = i — 1 and i — x are handled similarly, with the 
only difference that the rates include also the constant terms e^ or e~*. D 

Fix a path z{t) in Z with z(0) = and since t is fixed abbreviate z — z(t). 
Define 

r+l z(t) 

n=z(t) + l n=£-l 

Due to the frozen columns at the boundaries and the normalization /io(0) — 
we have the identity 

h,{t) = I+{t)-I+{0) = -/_(<) +/_(0) 

from which 

(A.9) Var[/i,(t)] = Cov[-/_(t) + /_(0) , I+{t) - /+(0)] . 

For each time t the increment variables {LOi{t) : £ < i < x} are i.i.d. /i^- 
distributed. The process is independent of the initial values cjf_i(0) and cjr+i(0) 
of the boundary increments. This is because the bricklayer at site i ~ 1 lays 
bricks to his right at rate e* regardless of the value of the increment at site £—1, 
and similarly the left action of the bricklayer at site r+l has constant rate e~*. 
But the time increment 

(A.IO) W£_i(t) - w£-i(0) = -he-i{t) + he-i{0) 
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and the same for io^^i(t) are needed for hz{t). Expanding and removing van- 
ishing covariances from (|A.9p leaves 

(A.ll) Var[/i,(i)]= Yl Cov(w,(0),w,(i)) + ^ Cov(cj,(i),c^,(0)) 
e<i<o e<i<z 

z<j<z 0<j<r 

(A.12) + Y. Cov(cj,;(0),w,+i(t)) + Y. Cov(w,_i(t),w,(0)) 

f<!<0 0<j<r 

(A.13) - Y Cov(w,(t),Lj,+i(i)) - ^ Cov(tj^_i(i),w,(t)) 

l<i<z z<j<t 

(A.14) - Cov(cj,_i(t),L^,+i(t)). 

We argue that the sums on lines (jA.12[) ~ (IA.14l) vanish as —ijV — > oo by 
showing that the covariances decay exponentially in the spatial distance. 

We illustrate with a term Cov(w.JQ(t),ajr+i(i)) for a fixed io ^ z from the 
first sum on line (IA.13[) . Consider — ^,r large so that £ < io < z < x. Let 
w = [(io + '^)/2j be a lattice point at or next to the midpoint between ip and 
r. Let ({t) and ^(i) be two further processes whose initial configurations agree 
with those of w(t), both for the increments and for the heights: 

(A.15) C(0) = C(0) = w(0) and h^{0) = h^{0) = h{0). 

Processes ({t) and ^(i) follow the {i, r, 9) evolution, except that the columns 
itkiit) : i > w ~ 1} with bases in [w — 1, oo) are not permitted to grow, and the 
columns {h^{t) : i < w} based in (—00,1/7] are similarly frozen. (Equivalently, 
replace the Poisson clocks of the corresponding bricklayers with empty point 
measures.) To determine their dynamics, in addition to disjoint collections of 
Poisson clocks, C{t) requires initial increments {cL'i(O) : i < w — 1} while ^{t) 
requires initial increments {^^(0) : i > w + 1}. Thus {Ci(i) '■ i ^ u; — 1} and 
{(,i{t) : i > w + 1} are independent processes. The height processes {h'^{t) : 
i < w — 2} and {h^{t) : i > w + 1} are not independent. For example if z > 
the initial heights {/if (0) : i > w + 1} are very much influenced by the initial 
increments {Ci(0) : < i < z}. 

On the sites where C{t) and £,{t) lay bricks they obey the same realizations 
of Poisson clocks as the (i, r, 6) process ui{t). This coupling preserves the in- 
equalities 

(A. 16) h^{t) < h{t) and h^it) < h{t). 

Let 

(A.17) A^{uj,{t)^Cr{t)ioTi<io} and B ^ K+i(t) = Cr+i(t)}- 

Let us say that the space-time window [a, 6] x (sq, si] is a block if hi(si) = 
/li(so) for some \a\ < i < [6J — 1, in other words, if some column failed to 
grow inside the space interval [a, b] during time interval {sq, si\. The space-time 
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window is not a block if every column inside grew during (sq, si]. The sense 
of the terminology is that a block does not permit a discrepancy to pass. For 
convenience we do not require a, b to be integers. We restate Corollary 5.5 from 
[^. This lemma is valid in the stationary process because there is control over 
the spatial averages of the increments uJi and thereby control over rates. 

Lemma A. 6. [S', Cor. 5.5] There exist constants /cq < oo and t^ > and a 

function < G{s) < oo for s £ (0,to) such that G{s) /* cx) as s \ 0, and this 
bound holds: if si ~ sq < to and b — a > ko then 

P{[a,b] X {sQ,si] is not a block} < e-(^-'^)-G('5i-so)_ 
For a given 9 this bound is valid for all {£, r, 9) processes. 

Recalling that t is fixed in the present discussion, fix a positive integer m 
and real r > so that 

t — mr and t e (0, to]. 

On the event A"^ there must exist a sequence of times < tw-2 < tw-3 < 
tw-4: < ■ ■ ■ < tig < t such that column hi grew at time ti, io < i < w — 2. This 
results from the observation that in the range i < w — 2 the first discrepancy hi — 
h^>Om column heights at i can appear only next to an existing discrepancy. 
Thus the leftmost discrepancy Q{t) — min{i : hi{t) — h'^{t) > 0} starts from 
the value (5(0) = oo due to (|A.15p . arrives at the boundary w — 1 of the frozen 
region of h at time twi-i when h^^^i first grows, and then moves to the left 
one step at a time and always with a jump of the column hqu^iit) that is not 
matched by a jump in hQ,^.^{t). 

Consequently event A'^ implies that at least one of the windows 

w — 2 — if) , ,1/7 — 2 — in 1 

-2-j °-,w~2-{j-l) "- x[{j-1)t,jt], l<j<m, 

m m J 

is not a block. For if all these windows were blocks, the sequence of growths 
over intervals [w ~ i^w — i -\- 1] at times im-i, 2 < i < w — ig, could not happen. 
This gives the bound 



PiA") < me"^(™"2-»o)/™ < 



me 



-C(t-io)/m 



where for the second inequality we assumed that t is large enough relative to z. 
The same argument for the propagation of discrepancies can be repeated for 
event B in (jA.17[) to improve the bound to 

(A.18) P(yl" U B") < 2me"^('^~'°)/™ < Cie~^^^' 



-io) 



where we subsumed m in the constants. 

Next we turn these bounds into bounds on covariances. 

Cov(w,„(i),c^,+i(i)) = Cov(C»„(i),et+i(i)) 

+ E[(w,„(tK+i(t) - C.o(i)6+iW)lAcus^ 

- E[(w,„(t) - Go W)1ac] • Ew,+i(i) 

- E[(a;,+i(t)-C,+i(t))lBe].EC.o(i). 
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By independence oi Qg{t) and ^r+i(0 the first covariance after tfie equality sign 
vanislies. We claim that for £ < i < x:, 

(A.20) E\io,{t)\P + E|C.(t)r + E|c^,+i(i)|P + E|6+i(i)r < e^*. 

Granting this for the moment, we apply Holder's inequality to the other expec- 
tations in (jA.igp . then (|A.18|) and the moment bounds (IA.20p to arrive at 

(A.21) \Cov{uj,,{t),u;,+i{t))\ < Cie^^^'^'-'^l 

Now to verify (|A.20p term by term. For Wi(t) this is simply a finite moment 
under the jjl^ distribution. The other increment variables we express so that 
Lcinma lA.5l applies to each case. For the boundary increment u]^^i{t) write 

cj,+i(t) = h,{t) - h,+i[t) = h,{t) - h,{Q) + h,{Q) - h,+i{Q) 
■ =/i,(f)-/i,(0)+a;,(0) 

utilizing the frozen column h-^+i as in (JA.IOI) . Utilizing (IA.15|) and (|A.16[) 

Ut) = hi_^{t) - h\{t) < /i,_i(t) - h,{0) = h,^i{t) - /i,_i(0) +cj,(0) 
with a similar upper bound for — Ci(i), which together give 

\Q{t)\ < {h,-i[t) - /i,_i(0)) V {h,{t) - h,{0)) + |u;,(0)|. 
For ^r-|_i(t) use stochastic monotonicity of columns to begin with 

e.+i(t) = hi[t) - hi^,{t) < Kit) - /i,+i(0) 

and continue as in (|A.22p . The completes the verification of (jA.20p . 

We have proved the exponential decay of the covariance in (|A.2ip . The 
same arguments work for all the covariances on lines (jA.lip - (jA.14p . and we 
state the lemma in this generality. The bounds extend also to the infinite- 
volume stationary process because Lemma (jA.Sp gives moment bounds that 
ensure uniform integrability. 

Lemma A. 7. There exist constants Ci = Ci(t,9) e (0, oo) such that, for all 
stationary (£, r, 9) processes and all i^jGZ and s £ [0,t], 

(A.23) Cov{uji{s),ujj{t)) < Cie-^=l*-^l. 

The bound is also valid for the infinite-volume stationary process. 

Now we can complete the proof of the covariance formula. 

Proof of Theorem \A.l[ Since the bound (JA.23P hold uniformly as —£,x — > oo, 
the sums on lines (jA.12p - (IA.14p vanish while the sums on line (jA.lip converge 
to 

^ Cov(^,(0),w,(i)) + Y. Cov(a;,(i),a;,(0)). 

i<0.,j>z i<z,j>0 

Translation invariance of the stationary infinite volume process turns the above 
sum into the right-hand side of (jA.ip . The left-hand side of (jA.lip converges 
to the left-hand side of (jA.ip as —£,1 — > oo by Lemma [A. 21 and by the uniform 
integrability given by Lemma IA.5I D 
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